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^L(- Abstract. The present paper proves that all limit points of sequences 

' of renormalized solutions of the Boltzmann equation in the limit of small, 

asymptotically equivalent Mach and Knudsen numbers are governed by 
Leray solutions of the Navier-Stokes equations. This convergence result 
holds for hard cutoff potentials in the sense of H. Grad, and therefore 
completes earlier results by the same authors [Invent. Math. 155, 81-161 
P_l ' (2004)] for Maxwell molecules. 
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1. Introduction 



The subject matter of this article is the derivation of the Navier-Stokes 
equations for incompressible fluids from the Boltzmann equation, which is 
^ ' the governing equation in the kinetic theory of rarefied, monatomic gases. 

On I In the kinetic theory of gases founded by Maxwell and Boltzmann, the 

' state of a monatomic gas is described by the molecular number density in 

. the single-body phase space, / = f{t,x,v) > that is the density with 

I respect to the Lebesgue measure dxdv of molecules with velocity v G R'^ 

• and position x £ R^ at time t > 0. Henceforth, we restrict our attention to 

CO , the case where the gas fills the Euclidian space R^. For a perfect gas, the 

' number density / satisfies the Boltzmann equation 



(1.1) dtf + v-V,f = B{f,f), x,vGTi\ 

where B{f, f) is the Boltzmann collision integral. 
5^ I The Boltzmann collision integral acts only on the v variable in the number 

density /. In other words, ;S is a bilinear operator defined on functions of 
the single variable v, and it is understood that the notation 

(1.2) B{f, f)it, X, v) designates ^(/(i, x, •), fit, x, ■)){v) , 

For each continuous / = f{v) rapidly decaying at infinity, the collision 
integral is given by 

(1.3) B{f,f)iv)= [[ {f{v')f{v[)-f{v)f{vi))biv-vi,u;)dvicLj 
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where 

v' = v' (v,Vi,Lil) = V — (v — Vl) ■ UJUJ , 

(14) V > . ; \ ; ^ 

Vi = Vi[V, Vl,UJ) = Vl + [V — Vl) ■ UJUJ . 

The colhsion integral is then extended by continuity to wider classes of 
densities /, depending on the specifics of the function b. 

The function b = b{v — vi,uj), called the collision kernel, is measurable, 
a.e. positive, and satisfies the symmetry 

(1.5) b{v — = b{vi — v,Lu) = b{v' — v[,uj) a.e. in {v,vi,ijj) . 

Throughout the present paper, we assume that b satisfies 

< b{z,uj) < Cb{l + \z\ f\ cos(f7w)| a.e. on x 

(1-6) f I \z\ o 

/ b(z,uj)doj > — -— a.e. on R . 

Js^ - Cbl + \z\ 

for some C;, > and /? € [0, 1]. The bounds (|1.6p are verified by all collision 
kernels coming from a repulsive, binary intermolecular potential of the form 
U{r) = Uo/r^ with Grad's angular cutoff (see fl5]) and s > 4. Such power- 
law potentials are said to be "hard" if s > 4 and "soft" otherwise: in 
other words, we shall be dealing with hard cutoff potentials. The case of a 
hard-sphere interaction (binary elastic collisions between spherical particles) 
corresponds with 

(1.7) b{z, LU) = \z ■ L0\ ] 

it is a limiting case of hard potentials that obviously satisfies p.6p . even 
without Grad's cutoff. At the time of this writing, the Boltzmann equation 
has been derived from molecular dynamics — i.e. Newton's equations of 
classical mechanics applied to a large number of spherical particles — in 
the case of hard sphere collisions, by O.E. Lanford [16], see also [9] for the 
case of compactly supported potentials. Thus the collision kernel b given by 
()1.7p plays an important role in the mathematical theory of the Boltzmann 
equation. 

The only nonnegative, measurable number densities / such that B{f, f) = 
are Maxwellian densities, i.e. densities of the form 

(1.8) f(v) = ^^-A_e^ =: MR,u,e{v) 

for some i? > 0, > and U G R^. Maxwellian densities whose parameters 
R, U, are constants are called "uniform Maxwellians" , whereas Maxwellian 
densities whose parameters R, U, Q are functions of t and x are referred to as 
"local Maxwellians". Uniform Maxwellians are solutions of (jl.ip : however, 
local Maxwellians are not solutions of (jl.ip in general. 

The incompressible Navier-Stokes limit of the Boltzmann equation can be 
stated as follows. 
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Navier-Stokes Limit of the Boltzmann Equation 

Let u*" = G be a divergence- free vector field on R'^. For eacli 

e > 0, consider the initial number density 

(1-9) fe^'ix^v) = • 

Notice that the number density /*" is a slowly varying perturbation of order 
e of the uniform Maxwellian A^i.o,i- Let solve the Boltzmann equation 
(jl.ip with initial data ()1.9p . and define 

(1.10) Ue{t,x):=- vfe ( ^,-,v] dv . 

Then, in the limit as e — > 0"*" (and possibly after extracting a converging 
subsequence), the velocity field satisfies 

Ue^u in P'(R+ X R^) 

where u is a solution of the incompressible Navier-Stokes equations 

dtu + divx{u (^u) + VxP = i^AxU , x e R^ , t > , 
dwx u = , 

with initial data 

(1.12) «lt=o = ^'"- 

The viscosity u is defined in terms of the collision kernel b, by some implicit 
formula, that will be given below. 

(More general initial data than (11. 9[) can actually be handled with our 
method: see below for a precise statement of the Navier-Stokes limit theo- 
rem.) 

Hydrodynamic limits of the Boltzmann equation leading to incompress- 
ible fluid equations have been extensively studied by many authors. See 
in particular [2\ for formal computations, and pj, |3j for a general program 
of deriving global solutions of incompressible fluid models from global so- 
lutions of the Boltzmann equation. The derivation of global weak (Leray) 
solutions of the Navier-Stokes equations from global weak (renormalized a 
la DiPerna-Lions) solutions of the Boltzmann equation is presented in [3], 
under additional assumptions on the Boltzmann solutions which remained 
unverified. In a series of later publications [201 [22l HI [10] some of these 
assumptions have been removed, except one that involved controlling the 
build-up of particles with large kinetic energy, and possible concentrations 
in the x-variable. This last assumption was removed by the second author 
in the case of the model BGK equation [231 [24], by a kind of dispersion 
argument based on the fact that relaxation to local equilibrium improves 
the regularity in v of number density fluctuations. Finally, a complete proof 
of the Navier-Stokes limit of the Boltzmann equation was proposed in [13]. 
In this paper, the regularization in v was obtained by a rather different 
argument — specifically, by the smoothing properties of the gain part of 
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Boltzmann's collision integral — since not much is known about relaxation 
to local equilibrium for weak solutions of the Boltzmann equation. 

While the results above holds for global solutions of the Boltzmann equa- 
tion without restriction on the size (or symmetries) of its initial data, earlier 
results had been obtained in the regime of smooth solutions [71 H]. Since 
the regularity of Leray solutions of the Navier-Stokes equations in 3 space 
dimensions is not known at the time of this writing, such results are limited 
to either local (in time) solutions, or to solutions with initial data that are 
small in some appropriate norm. 

The present paper extends the result of [13] to the case of hard cutoff 
potentials in the sense of Grad — i.e. assuming that the collision kernel 
satisfies ()1.6p . Indeed, [13] only treated the case of Maxwell molecules, for 
which the collision kernel is of the form 



The method used in the present paper also significantly simplifies the original 
proof in [13] in the case of Maxwell molecules. 

Independently, CD. Levermore and N. Masmoudi have extended the anal- 
ysis of fL3\ to a wider class of collision kernels that includes soft potentials 
with a weak angular cutoff in the sense of DiPerna-Lions: see [T7]. Their 
proof is written in the case where the spatial domain is the 3-torus R^/Z^. 

In the present paper, we handle the case of the Euclidian space R^, which 
involves additional technical difficulties concerning truncations at infinity 
and the Leray projection on divergence-free vector fields — see Appendix C 
below. 



2.1. Global solutions of the Boltzmann equation. The only global ex- 
istence theory for the Boltzmann equation without extra smallness assump- 
tion on the size of the initial data known to this date is the R. DiPerna-P.-L. 
Lions theory of renormalized solutions [HJ [18] . We shall present their the- 
ory in the setting best adapted to the hydrodynamic limit considered in the 
present paper. 

All incompressible hydrodynamic limits of the Boltzmann equation in- 
volve some background, uniform Maxwellian equilibrium state — whose role 
from a physical viewpoint is to set the scale of the speed of sound. With- 
out loss of generality, we assume this uniform equilibrium state to be the 
centered, reduced Gaussian density 



b{z, LO) 



cos(2:,cj)|6*(| cos{z,uj)\) with — < 6=,, < C=k . 



2. Formulation of the problem and main results 



(2.1) 



M{v) ■.= Mi,o,i{v) 



1 




(27r)3/2 
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Our statement of the Navier-Stokes limit of the Boltzmann equation given 
above suggests that one has to handle the scaled number density 

(2.2) F,{t,x,v) = fj^,-,v 



where /g is a solution of the Boltzmann equation (jl.ip . This scaled number 
density is a solution of the scaled Boltzmann equation 

(2.3) e'^dtF, + ev ■ V^F, = B{F„ F,) , x, t; G R^ t > . 

Throughout the present section, e is any fixed, positive number. 

Definition 2.1. A renormalized solution of the scaled Boltzmann equation 
112. 3\) relatively to the global equilibrium M is a function 

FgC(R+,L1„,(R3xR3)) 

such that 

r' (^) B{F, F) G lL(R+ X r3 X r3) 
and which satisfies 

(2.4) M {e'd, + ev • V.) T (^) = T' (^) B{F, F) 

for each normalizing nonlinearity 

r G C\'R+) such that \T'{z)\ < ^ ^ , z>0. 

Vl + 

The DiPerna-Lions theory is based on the only a priori estimates that 
have natural physical interpretation. In particular, the distance between any 
number density F = F{x, v) and the uniform equilibrium M is measured in 
terms of the relative entropy 

(2.5) H{F\M):= jj \n (^^^ - F + dxdv . 
Introducing 

(2.6) h{z) = {l + z)\n{l + z)-z>Q, z > -I , 
we see that 

H(F\M) =11 h (— - 1 ) Mdvdx > 

with equality if and only ii F = M a.e. in x, v. 

While the relative entropy measures the distance of a number density F to 
the particular equilibrium M, the local entropy production rate "measures 
the distance" of F to the set of all Maxwellian densities. Its expression is as 
follows: 

(2.7) S{F) = 1 [[[ (F'Fi'-FFi)ln ( ^) b{v-v,,u;)dvdviAj . 
The DiPerna-Lions existence theorem is the following statement [H [18] . 
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Theorem 2.2. Assume that the collision kernel b satisfies Grad's cutoff 
assumption ^.6\) for some (3 £ [0, 1]. Let F*" = F*"(x, v) be any measurable, 
a.e. nonnegative function on x R'^ such that 

(2.8) /7(F™|Af) < +CX). 

Then, for each e > 0, there exists a renormalized solution 

eC7(R+,lL(RxR3)) 

relatively to M of the scaled Boltzmann equation i2.3\) such that 

p I — 



Moreover, satisfies 

(a) the continuity equation 

(2.9) edt I F.du + div^ / 



vF^dv 

R3 JR3 



and 

(b) the entropy inequality 

ft 



(2.10) H(FJM)(t) + \ f [ £(F,)(s,x)dsdx<H(FnM), t>0. 

Jo Jr3 



Besides the continuity equation (j2.9p . classical solutions of the scaled 
Boltzmann equation (j2.3p with fast enough decay as — > oo would satisfy 
the local conservation of momentum 

(2.11) edt vF^dv + div^ v0vF^dv = O, 

JR3 JR3 

as well as the local conservation of energy 

(2.12) edt ^\vfF,dv + dw^ v^\v\'^F,dv = . 

JR3 yR3 



Renormalized solutions of the Boltzmann equation (j2.3p are not known to 
satisfy any of these conservation laws except that of mass — i.e. the continu- 
ity equation (j2.9p . Since these local conservation laws are the fundamental 
objects in every fluid theory, we expect to recover them somehow in the 
hydrodynamic limit e ^ 0"^. 

2.2. The convergence theorem. It will be more convenient to replace the 
number density F^ by its ratio to the uniform Maxwellian equlibrium M; 
also we shall be dealing mostly with perturbations of order e of the uniform 
Maxwellian state M. Thus we define 

F G - 1 

(2.13) G, = 9, g, = ^^^. 

Likewise, the Lebesgue measure dv will be replaced with the unit measure 
Mdv, and we shall systematically use the notation 

(2.14) (</>)=/ (l){v)M{v)dv , for each cl) e (Mdv) . 
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For the same reason, quantities like the local entropy production rate involve 

the measure 

(2.15) 



dfi{v,vi,uj) = b{v — vi,uj)MidviMdvduj , /// dfi{v,vi,u;) = 1 , 

whose normalization can be assumed without loss of generality, by some 
appropriate choice of physical units for the collision kernel b. We shall also 
use the notation 
(2.16) 

(ip) = /// 'tp{v,vi,uj)dn{v,vi,uj) for G L^(R^ xR"^ X S'^,dn) . 

i J JR3xR3xS2 

From now on, we consider solutions of the scaled Boltzmann equation 
()2.3p that are perturbations of order e about the uniform Maxwellian M. 
This is conveniently expressed in terms of the relative entropy. 

Proposition 2.3 (Uniform a priori estimates). Let F^*" = F^^{x,v) be a 
family of measurable, a.e. nonnegative functions such that 

(2.17) sup \h{F^''\M) = C^" < +00 . 

Consider a family (F^) of renormalized solutions of the scaled Boltzmann 
equation \2. 3)) with initial data 

(2.18) 

Then 

(a) the family of relative number density fluctuations satisfies 

(2.19) 1 I {h{eg,{t,x,-)))dx<C'^ 

JR3 

where h is the function defined in \2. 0(1 : 

(b) the family ^{y/Gl — 1) is bounded in L'^iYl^; L'^{Mdvdx)): 

(2.20) / / f^/^I^y\d^<c" ; 



u=o 



/R3 \ V e 

(c) hence the family g^ is relatively compact in Lj^^{dtdx; {Mdv)) ; 

(d) the family of relative number densities satisfies the entropy pro- 
duction — or dissipation estimate 

Proof. The entropy inequality implies that 



H{F,\M){t)= [ {h{G,-l)){t,x)dx <H{Fi''\M) <&''e^ 

JR3 



'R3 ' 

which is the estimate (a). 
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The estimate (b) follows from (a) and the elementary identity 

h{z-l) - i^-lf = zlnz- - (^/I- 1)2 

= 2zln^/I-2{^/^-l)^/I 

= 2V^(V^lnV^- 1) > 

From the identity 



(2.22) g, = 2^—^ + e 

and the bound (b), we deduce the weak compactness statement (c). 
Finally, the entropy inequality implies that 



/ / £{F,){s,x)dxds <&''e^ . 
Jo Jr3 



£{F,) = 1111 {F'Xi - FeF.i) In ( ] b{v - vuu;)dvdvidio 



GeGtl 



Observing that 

kf: 

/R3xR3xS2 \FeF^i 

= i((G',G',i-aG'ei)ln 
and using the elementary inequality 

l(X-Y)lny>iVX-VYf, X,Y>0 
leads to the dissipation estimate (d). □ 



Our main result in the present paper is a description of all limit points of 
the family of number density fluctuations g^. 

Theorem 2.4. Let Fl^ he a family of measurable, a.e. nonnegative func- 
tions defined on R'^ x R'^ satisfying the scaling condition {2.11^. Let he 



a family of renormalized solutions relative to M of the scaled Boltzmann 
equation \2. 3|) with initial data 12.18\) . for a hard cutoff collision kernel b 
that satisfies il.6\} with (3 G [0,1]. Define the relative numher density 
and the numher density fluctuation g^ hy the formulas i2.13\) . 

Then, any limit point g in L\^^{dtdx; L^ [M dv)) of the family of numher 
density fluctuations g^ is an infinitesimal Maxwellian of the form 

g{t, X, v) = u{t, x) ■ V + 6{t^x)\{\v\^ — 5) , 

where the vector field u and the function 9 are solutions of the Navier-Stokes- 
Fourier system 

dtu + diva;(n ®u)+ V xP = v^xU , divx- u = 

^^■^^^ dtO + A\wx{ue) = kI:^xO 
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u; - lim P ( - / vFTdv 



with initial data 



(2.24) 

= u; - lim ^ j {\\v\'^ - 1)(F,*" - M)dv , 

where P is the Leray orthogonal projection in L^(R^) on the space of diver- 
gence-free vector fields and the weak limits above are taken along converging 
subsequences. Finally, the weak solution {u, 9) of ^2.23\) so obtained satisfies 
the energy inequality 

ft 



(2.25) 



( {\\u{t,x)\'^ + \\9{t,x)\'^)dx+ [ [ {iy\V^u\'^ + ^K\V^e\'^)dx 

R3 Jo Jrs 



< lim \h(F^''\M) 
e->o+ e 



The viscosity v and thermal conductivity k are defined implicitly in terms of 
the collision kernel b by the formulas 1^2.27 ) below. 



There are several ways of stating the formulas giving i' and k. Perhaps 
the quickest route to arrive at these formulas is as follows. 

Consider the Dirichlet form associated to the Boltzmann collision integral 
linearized at the uniform equilibrium M: 

(2.26) Vm{^) := 1(1^' + . 

The notation | • p designates the Euclidian norm on R'^ when $ is vector- 
valued, or the Frobenius norm on M^CR) (defined by \A\ = tiace{A* A)^^'^) 
when $ is matrix- valued. Let D* be the Legendre dual of P, defined by the 
formula 

P*(^') :=sup((^'-$) -P($)) 

where the notation ^(v) ■ '^(v) designates the Euclidian inner product in 
whenever ^ are vector valued, or the Frobenius inner product in Af3(R) 
whenever $, ^ are matrix- valued (the Frobenius inner product being defined 
hy A - B = tmce{A*B).) 

With these notations, one has 

(2.27) r,:=lv*{v(S)v-l\v\^I), n := ^V*{^v{\v\^ - 5)) . 

The weak solutions of the Navier-Stokes-Fourier system obtained in The- 
orem [23] satisfy the energy inequality ()2.25p and thus are strikingly similar 
to Leray solutions of the Navier-Stokes equations in 3 space dimensions — 
of which they are a generalization. The reader is invited to check that, 
whenever the initial data F^*" is chosen so that 

^ii-(F™|M) ^ i / \u'''{x)\'^dx as e ^ 0+ , 
e Jr3 
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then the vector field u obtained in Theorem 12.41 is indeed a Leray solution 
of the Navier-Stokes equations. More information on this kind of issues 
can be found in [13]. See in particular the statements of Corollary 1.8 and 
Theorem 1.9 in [13], which hold verbatim in the case of hard cutoff potentials 
considered in the present paper, and which are deduced from Theorem 12.41 
as explained in |13j . 



2.3. Mathematical tools and notations for the hydro dynamic limit. 

An important feature of the Boltzmann collision integral is the following 
symmetry relations (the collision symmetries). These collision symmetries 
are straightforward, but fundamental consequences of the identities (jl.Sp 
verified by the collision kernel, and can be formulated in the following man- 
ner. Let $ = $(v,ui) be such that G L^iYi? x x S^,d^). Then 
(2.28) 




^{v,vi)diJL{v,vi,u:) = III ^{vi,v)dfi{v,vi,uj) 
R3xR3xS2 JJJr3xR3xS2 

^{v'{v, vi,u;),v[{v, vi,Lu))dn{v, Vi,tj) 

R3xR3xS2 

where v' and v'l are defined in terms of v,vi,uj by the formulas (jl.4p . 

Since the Navier-Stokes limit of the Boltzmann equation is a statement 
on number density fluctuations about the uniform Maxwellian M, it is fairly 
natural to consider the linearization at M of the collision integral. 

First, the quadratic collision integral is polarized into a symmetric bilinear 
operator, by the formula 

BiF, G) := i {B{F + G,F + G)- B{F, F) - B{G, G)) . 

The linearized collision integral is defined as 

(2.29) Cf = -2M~^B{M,Mf) . 

Assuming that the collision kernel b comes from a hard cutoff potential 
in the sense of Grad (jl.6p . one can show (see [15] for instance) that C is 
a possibly unbounded, self-adjoint, nonnegative Fredholm operator on the 
Hilbert space L'^(R^ , Mdv) with domain 

D{£) = L'^(R^,a{\v\fMdv) 

and nullspace 

(2.30) Ker £ = span{l, fi, f2, fs, , 
and that C can be decomposed as 

£g{v) = a{\v\)g{v) - lCg{v) 

where /C is a compact integral operator on Lp'{Mdv) and a = a{\v\) is a 
scalar function called the collision frequency that satisfies, for some C > 1, 

0<a_ <a{\v\) <a+{l + \v\f . 
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In particular, L has a spectral gap, meaning that there exists C > such 
that 



{Madv) 



(2.31) (/£/)> c||/-n/||i, 

for each / E D[C), where 11 is the orthogonal projection on Keri2 in 
L2(R3,Mdf), i.e. 

(2.32) n/ = (/) + {vf) . V + (dlH' - l)f)\{\v\^ - 3) . 

The bilinear collision integral intertwined with the multiplication by M 
is defined by 

(2.33) Q{f,g)=M-^B{Mf,Mg). 

Under the only assumption that the collision kernel satisfies (jl.Sp together 
with the bound 

(2.34) / b{z,uj)duj < a+{l + \z\)f^ , 

Q maps continuously L'^(R?,M{1 + \v\)f^dv) into L'^ (Ti? , a^^ M dv) . Indeed, 
by using the Cauchy-Schwarz inequality and the collision symmetries (j2.28p 
entailed by (11.5(1 
(2.35) 

ll2(5,Mlli.(.-iMd.) = ^3«(lH)-' 

X J J {g'h[ + g[h' - ghi - gih)h{v - vi,uj)Midviduj^ Mdv 
< \ I a{\v\y^ [ II b{v -vi,uj)Midviduj 



R3 

{g'h\ + g'lh' — ghi — gih)^b{v — vi,uj)AIidviduj^ Mdv 
< sup a(|u|)^^ // b{v — vi,uj)Midviduj 

{{g'h[f + {g[h'f + igh.f + ig,hf)d^i{v,vi,u;) 




2 I /„ u\2\ 



<2C {{ghiY + {gihy) { b{v -vi,uj)duj] MMidvdvi 

Another important property of the bilinear operator Q is the following re- 
lation: 

(2.36) Q(/, /) = for each / G Ker £ , 

which follows from differentiating twice both sides of the equality 
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with respect to i? > 0, 6 > and [/ G — see for instance flas (59-60) 
for a quick proof of this identity. 

Young's inequahty 

Since the family of number density fluctuations satisfies the uniform 
bound (a) in Proposition 12.31 and the measure Mdv has total mass 1, the 
fluctuation can be integrated against functions of v with at most quadratic 
growth at infinity, by an argument analogous to the Holder inequality. This 
argument will be used in various places in the proof, and we present it here 
for the reader's convenience. To the function h in (j2.6p . we associate its 
Legendre dual h* defined by 

h*(C) ■= sup (C^ - h{z)) = - C - 1 . 

Thus, for each C > ^^'^ each z > —1, one has 

(2.37) C\z\ < h{\z\) + h*{C) < h{z) + h*{C) 

since 

h{\z\)<h{z), z>-l. 

The inequality (I2.37P is referred to as the Young inequality (by analogy with 
the classical Young inequality 

zP 

Cz<— + — , z,C>0 
P Q 

which holds whenever 1 < p,q < oo satisfy ^ + ^ = 1.) 

Notations regarding functional spaces 
Finally, we shall systematically use the following notations. First, Lebesgue 
spaces without mention of the domain of integration always designate that 
Lebesgue space on the largest domain of integration on which the measure 
is defined. For instance 

I/iMdv) designates LP (R^Mdw) 
LP{Mdvdx) designates LP(R^ x R^; Mdvdx) 
LP{dfi) designates LP(R^ x R^ x S^; dfi) 

When the measure is the Lebesgue measure, we shall simply denote 
LP := LP{K^; dx) , Lf_^ := LP(R+ x R^; dtdx) . 

Whenever is a normed space, the notations 0{5)e and o{5)e designate 
a family of elements of E whose norms are 0{6) or o{5). (For instance 0{\)e 
designates a bounded family in while 0(1)^; designates a sequence that 
converges to in E.) 

Although L^^^ spaces are not normed spaces, we designate by the notation 
0[5)iv (f^) a family G L^^JJl) such that, for each compact K <Z 

WfeUnK) = 0{5) . 
The notation o{6)tp /q-, is defined similarly. 

loc ^ ' 



NAVIER-STOKES LIMIT OF THE BOLTZMANN EQUATION 13 

2.4. Outline of the proof of Theorem 12.41 In terms of the fluctuation 
Qe, the scaled Boltzmann equation (j2.3p with initial condition ()2.18p can be 
put in the form 

edtge + V ■ V^g, = --C{ge) + Q{ge,ge) , 

ge\t=o = gT ■ 

Step 1 : We first prove that any limit point g of the family of fluctuations 
as e — > 0"*" satisfies 

g = ^g 

where 11 is the orthogonal projection on the nullspace of C defined in (j2.32p . 

Hence, the limiting fluctuation g is an inflnitesimal Maxwellian, i.e. of 
the form 

(2.39) g{t, X, v) = p{t, x) + u{t, x) ■ v + 9{t, x)^{\vf - 3) . 

The limiting form of the continuity equation (j2.9p is equivalent to the in- 
compressibility condition on u: 

divx u = . 

Step 2 : In order to obtain equations for the moments 

P={g), u = {vg), and 9 = {{^\v\'^ - l)g) 

we pass to the limit in approximate local conservation laws deduced from 
the Boltzmann equation in the following manner. 

Besides the square-root renormalization, we use a renormalization of the 
scaled Boltzmann equation (|2.3|) based on a smooth truncation 7 such that 

(2.40) 7GC-(R+,[0,1]), 7|[oJ, = l, 7|[2,+oo)=0- 
Define 

(2.41) = 
Notice that 

(2.42) supp(7) C [0,2], 7I ^ 3, = 1, and ||1-7||l- < 1 + IIt'IIl- • 

We use below the notation 7^ and 7^ to denote respectively 7(Ge) and 7(Ge). 

We also use a truncation of high velocities, defined as follows: given k > 6, 
we set 

(2.43) = k|lne|. 

For each continuous scalar function, or vector- or tensor-field ^ = ^(w), we 
denote by S^k^ the following truncation of S,: 

(2.44) ^KAv)=^ivn\v\'<K,- 
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Renormalizing the scaled Boltzmann equation ()2.3p with the nonlinearity 
T{Z) = (Z — l)'y{Z), we arrive at the following form of (j2.38p 

dt{ge%) + -V ■ V,(<7e7e) = ^%Q(.G,, G,) . 

Multiplying each side of the equation above by S^k^ , and averaging in the 
variable v leads to 

(2.45) dti^K^g^je) + div, -Jv^K^geje) = ^{CkMG'.G',, - G,G,i)} . 

Henceforth we use the following notations for the fluxes of momentum or 
energy: 

(2.46) F,(C) = -^{CK,9ele) 
with 

Civ) = A{v) := v®^ - \\v\^I , or C{v) = B{v) := \v{\v\^ - 5) . 
Likewise, we use the notation 

(2.47) D,(C) = ^(iKMG'fi',^ - GeGel)) 

for the conservation defect corresponding with the (truncated) quantity ^ = 
i{v), where ^ G span{l, vi,V2,v^, I^P}- 

The Navier-Stokes motion equation is obtained by passing to the limit as 
e ^ modulo gradient fields in the equation (j2.45p for ^{v) = Vj, j = 1, 2, 3, 
recast as 

(2.48) dt{vK,g,-f,) + div^ F,{A) + V^(i|t;||^^5e7e) = ^e{v) , 

while the temperature equation is obtained by passing to the limit in that 
same equation with ^{v) = ^(jvp — 5), i.e. in 

(2.49) dtililvl^ - f>)K.gele) + div, F,(S) = D,(i(|t;|2 - 5)) . 

For the mathematical study of that limiting process, the uniform a priori 
estimates obtained from the scaled entropy inequality are not sufficient. Our 
first task is therefore to improve these estimates using both 

a) the properties of the collision operator (see Section 3), namely a suitable 
control on the relaxation based on the coercivity estimate (j2.3ip 

{<PC<P)>CU-mhi^Mad.)-^ 

b) and the properties of the free transport operator (see Section 4), namely 
dispersion and velocity averaging. 

With the estimates obtained in Sections 3-4, we first prove (in Section 5) 
that the conservation defects vanish asymptotically 

De(C) ^ in L]^^{dtdx) , ^ G span{-ui, i;2, vs, \v?} 
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Next we analyze the asymptotic behavior of the flux terms. This requires 
splitting these flux terms into a convection and a diffusion part (Section 6) 

F,(C) - 2^C (n ^^~^ ^ ^ + |^CQ(yG:, ^ in LlMtdx) 

where C, is the unique solution in (Ker £)-*" of the Fredholm integral equation 

£C = C. 

For instance, the tensor field A and the vector field B defined by 

(2.50) A{v) := V ® V - \\v\'^ I , B{v) := \{\v\^ - f>)v 
satisfy 

(2.51) ^±Ker£, B_LKer£ 

componentwise, so that there exists a unique tensor field A and a unique 
vector field B such that 

(2.52) CA = A, CB = B, iandS_LKer£, 

The diffusion terms are easily proved to converge towards the dissipation 
terms in the Navier-Stokes-Fourier system 

2 

The formulas (j2.27p for the viscosity v and heat conduction k are easily 
shown to be equivalent to 

(2.53) j, = ^{A:A), ^ = ^{B.B). 

The (nonlinear) convection terms require a more careful treatment, involving 
in particular some spatial regularity argument and the filtering of acoustic 
waves (see Section 7). 

3. Controls on the velocity dependence 
of the number density fluctuations 

The goal of this section is to prove that the square number density fluc- 
tuation — or more precisely the following variant thereof 

2 



/a - 1 



CQ{^/G;, v^) ) ^ {C{v ■ V^g)) in Ll,{dtdx) . 



is uniformly integrable in v with the weight (1 + \ v\)p for each p < 2. 

In our previous work [13], we obtained this type of control for p = 
only, by a fairly technical argument (see Section 6 of [l3]). Basically, we 
used the entropy production bound to estimate some notion of distance 
between the number density and the gain part of a fictitious collision integral. 
The conclusion followed from earlier results by Grad and Caflisch on the v- 
regularity of the gain term in Boltzmann's collision integral linearized at 
some uniform Maxwellian state. 
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Unfortunately, this method seems to provide only estimates without the 
ht (l+|u|)^ (with/3 as in (jl.6p ) that is crucial for treating hard potentials 
other than the case of Maxwell molecules. Obtaining the weighted estimates 
requires some new ideas presented in this section. 

The first such idea is to use the spectral gap estimate (I2.31j) for the lin- 
earized collision integral. Instead of comparing the number density to (some 
variant of) the local Maxwellian equilibrium — as in the case of the BGK 
model equation, treated in \23i i24j . or in the case of the Boltzmann equa- 
tion with Maxwell molecules as in [13] — we directly compare the number 
density fluctuation to the infinitesimal Maxwellian that is its projection on 
hydrodynamic modes. 

The lemma below provides the basic argument for arriving at such esti- 
mates. 



Lemma 3.1. Under the assumptions of Theorem \2.4\ one has 
(3.1) 



/a 



1 



n 



L^iMdv) 



<0{e)L2^ + 0{e) 



/a 



1 



L^{Mdv) 



Proof. In order to simplify the presentation we first define some fictitious 
collision integrals C and Q 



R3xS2 



{g + gi- g' - g'i)Mih{v - vi,Lj)dvidw , 



Q{g.h) 



{g'h'i + g'lh' — ghi — gih)Mih{y — vi,uj)dvidu), 



'R3xS2 

obtained from C and Q by replacing the original collision kernel b with 

b{z, w) 



b{z, LO) 



1 + 



b{z, tJi)dLOi 



Start from the elementary formula 



(3.2) C 



/a 



1 



Q 



/a 



1 VGi-i 



Multiplying both sides of this equation by (I — n)(v^G7 — 1) and using the 
spectral gap estimate (I2.3ip leads to 



/a - 1 



n- 



1 



L'^(Mdv) 



(3.3) 



< e 



+ e 



/a 



Q 



1 ~ — — 



1 VGl-l 



L^{Mdv) 



L'^iMdv) 
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Denote 

dfjL{v,vi,uj) = MMib{v — vi,uj)dujdvdvi . 
By definition of b , one has 



b{v — vi,u>)du! < 1 . 
Hence Q is continuous on L'^{Mdv): by (j2.35p 

l|Q(5) ^)llL2(Mdi,) < '^\\9\\L'2(Mdv)\\HL'2(Mdv) ■ 

(Notice that b verifies (jl.Sp as does b). 
Plugging this estimate in (j3.3p leads to 



(3.4) 



/a 



1 



n- 



L^{Mdv) 



< Ce 





2 


+ e 


e 


L^{Mdv) 



L^{Mdv) 



Finally, applying the Cauchy-Schwarz inequality as in the proof of (|2.35p , 
one finds that 

2 



1 



L'^{Mdv) 



< sup 



eR3 J JR3xS2 



b{v — vi,u>)Midviduj 



< 



R3xR»xS2 



R3xR3xS2 



[V,Vl,UJ) 



since < 6 < 6. By the entropy production estimate (d) in Proposition 12.3 
the inequality above implies that 



L^iMdv) 



Oil] 



j2 . 



This estimate and p.4p entail the inequality p.l 



□ 



Notice that we could have used directly C and Q instead of their truncated 
analogues C and Q, obtaining bounds in weighted L? spaces by some loop 
argument, unfortunately much more technical than the proof above. 

The main result in this section — and one of the key new estimate in this 
paper is 



Proposition 3.2. Under the assumptions of Theorem \2.4\ for each T > 0, 
each compact K C R^, and each p < 2, the family 
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is uniformly integrable in v on [0, T] x K x R'^ with respect to the measure 
dtdxMdv. (This means that, for each 77 > 0, there exists a > such that, 
for each measurable ip = '^{x, v) verifying 

one has 



T 



^{i + \v\Y 



/a 



/O JKJH? 

Proof. Start from the decomposition 
(3.5) 



Mdvdxdt <r/.) 



J: = {l + \v\y 
+ (1 + \v\)^ 



1 



e 

/a - 1 



/a - 1 



[l + \v\)2 



{i + \v\yu 



n 



1 



e 

/a-1 



We recall from the entropy bound (b) in Proposition 12.31 that 



o(r 



L^{L'^{dxMdv)) 



SO that, by definition (j2.32p of the hydrodynamic projection 11 



(3.6) 



n 



Oil)Lf^(^(L'!{Mdv))) 



for all q < +00. Therefore the first term in the right-hand side of ()3.5 
satisfies 



(3.7) 



:i + \v\ 



n- 



1 



0{l)L^{Ll(Lr(A4dv))) 



for all < p < +00 and 1 < r < 2. 



In order to estimate the second term in the right-hand side of (|3.5p . we 
first remark that, for each 6 > 0, each p < 2 and each q < +00, there exists 
some C = C{p,q,5) such that 
(3.8) 



{i + \v\Yi^ 



/a-i 



(i + bl 



/a 



1 



0{5)Lf'{L'^{dxMdv)) +0 

donE 

:i + bi)p 



<:7(p,g,,5) 



L- (L9(M<ii,)) 



Indeed, by Young's inequality and Proposition 12.31 (a). 

2 r2 



1 



52 



52 /(l + lt;!? 



0(5^) 



L^iL^dxMdv)) 



52 

52 /(l + |7;|)f 
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We next use p.Sp with the two fohowing observations: first, the obvious 
continuity statement (13. 6p . Also, because of (j3.ip and the entropy bound 
(b) in Proposition 12.31 one has 



(3.9) 
Hence 



n- 



/a - 1 



L^{Mdv) 



0(6) 



L]^^(dtdx) 



{l + \v\f. 



{i + \v\y- 



<th{G,-lfl\l + \v\Yl^ 



/q- 1 

/a - 1 



n 
n 
n 



e 

- 1 



= : II + 111. 



Now 



II <^h{G,-l) + 6{l + \v\r 



n- 



+ 6{i + \v\y 



On the other hand 



\^ ^ mLj^^{dtdx;L^ (Mdv)) — ^ 



L^{L'^{Mdvdx)) + 0('^)L^(Li(AMi;da;)) 
(l + l^'lf 



+ 5 J 



(1 + |?;|)f/2exp 



2^2 



Li(Mdv) 



n 



/a -1 



Ll^^{dtdx;L2(Mdv)) 

= 0{5C{p,q,5)) 



with r = 

Putting ah these controls together shows that 



J </ + // + /// = 0(1) 



{Ll{L-{Mdv))) 



(3.10) 



+0((5)L^°o(^i(jVfdt,dx)) + 0{5)i,^(^ii(^Mdvdx)) + '^-'^ 

+0(5C(p,g,5))ii 

^{dtdx;V{Mdv)) 



I.e. 



(1-5)(1 + |H)^ 



/a - 1 



< 0(1)l^{L1(L'-(AMi)))) 



l^^{dtdx;L^{Mdv)) + 0(5)^^°° (j^^^^^^-)) 

which entails the uniform integrability in v stated in Proposition 



□ 
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Remark: replacing the estimate for II above with 

n<^h{G,-i) + ^-{i + \v\r 

and choosing 5 = 4 in (j3.1Up shows that 
/ /Q~ _ 1 \ ^ 

(1 + \v\f i j is bounded in L}^^{dtdx] L^{Mdv)) . 

In [3], the Navier-Stokes Umit of the Boltzmann equation is estabhshed as- 
suming the uniform integrabiUty in [0, T]xKx R'^ for the measure dtdxMdv 
of a quantity analogous to the one considered in this bound. As we shall 
see, the Navier-Stokes-Fourier limit of the Boltzmann equation is derived in 
the present paper by using only the weaker information in Proposition 13.21 

4. Compactness results for the number density fluctuations 
The following result is the main technical step in the present paper. 

Proposition 4.1. Under the assumptions in Theorem \2.4\ for each T > 0, 
each compact K C R^ and each p < 2, the family of functions 

is uniformly integrable on [0,T] x K x R^ for the measure dtdxMdv. 

This Proposition is based on the uniform integrability in v of that same 
quantity, established in Proposition 13.21 together with a bound on the 
streaming operator applied to (a variant of) the number density fluctua- 
tion (stated in Lemma 14. 2p . Except for some additional truncations, the 
basic principle of the proof is essentially the same as explained in Lemma 
3.6 of [T3] (which is recalled in Appendix B). In other words, while the result 
of Proposition 13.21 provides some kind of regularity in v only for the number 
density fluctuation, the bound on the free transport part of the Boltzmann 
equation gives the missing regularity (in the x-variable). 

The technical difficulty comes from the fact that the square-root renor- 
malization T{Z) = \fZ is not admissible for the Boltzmann equation due 
to the singularity at Z = 0. We will therefore use an approximation of the 
square-root, namely z \/ z + e°' for some a 2[. 

Lemma 4.2. Under the assumptions in Theorem \2. 41 for each a > 0, one 
has 

{edt + V ■ V,)^^^±^^ = 0(e2-/2)^,(^,„) 

+ 0{l)L2(^(^i_^_\^\)~i3Mdvdxdt) 

\^^{dtdx;L'^{{l+\v\)-l^ Mdv)) ■ 



n- 



+ -{i + \v\y 



/a 
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Proof of Lemma \4-S\ Start from the renormalized form of the scaled Boltz- 
mann equation (12. 3p . with normahzing function 



This equation can be written as 

(4.1) {edt + V . V.) ^'" ^^'-^ = l-^=i==Q(a, Ge) = Q\ + Ql 



with 

1 1 



' e2 2Ve" + G, 

2 



X // i^G'eG'a - Vaai) b{v - vuuj)dujMidvi 
(4.2) '''' ^ ^ 



X jj VGeGel (a/gJg^ - b{v - Vl,Vo)dujMidvi . 

The entropy production estimate (d) in Proposition 12.31 and the obvious 
inequahty 

Ve" + G, > e"/2 

imply that 

(4-3) WQlWi'^iMdvdxdt) < ^C*"e^ . 

On the other hand 



Write 
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Apply the Cauchy-Schwarz inequality as in the proof of ()2.35|) : then 



// 



^^Vg^g^_Vg^^^^ _ ^^^^^^^^^^^ 



L'^{{l+\v\)-(^Mdv) 



< sup((l + \v\)~'^l b{v - vi,uj)Midviduj 



1/2 



+ e sup + // Afi 



1/2 



h{v — vi,u!)dviduj 



1/2 



Therefore 



11 \/G^i^ — h{v - vi,uj)Midvi 



L'\(l+\v\)-I^Mdv) 



<C\ 1+e 



L2(Afi(l+|l,i|/3)dl,i)^ 



^ ^ei — \/ G^G 



el 



1/2 



because of the upper bound in Grad's cut-off assumption (jl.6p . 

Hence, on account of Proposition 13. 21 and the entropy production estimate 
(d) in Proposition! 



(4.4) Qi = 0(1)^2 

L\^^{dtdx;L'^ {{l+\v\)-P Mdv))) ■ 

Both estimates (|4.3p and ()4.4p together with ()4.ip entail the control in 
Lemma mH □ 



With Lemma 14.21 at our disposal, we next proceed to the 
Proof of Proposition \4-l\ 



Step 1. We claim that, for a > 1, 
/ ye" + - i ^ ^ 

(4.5) 



G, - 1 



— 0(e° ^)L^°°(L2^^((ia:;L2(Md»;))) 
+ 0(e"/^)ioo(ii(M(ij;(ia;)) • 
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Indeed, 



yje^ + G\ - 1 ^^G,-\ 



< 



G,>l/2 



+ e 1Ge<i/2 



(4.6) 



/a - 1 



and we conclude with the decomposition 



Ve" + a - 1 



/a - 1 



together with the fluctuation control (b) in Proposition 



/a - 1 



Step 2. Let 7 be a smooth truncation as in (I2.40|) . and set 



We claim that, for each fixed 6 > 0, 



(4.7) 
Indeed, 



{edt + v-V,)ct>i = 0{- ^ 



L]^^{Mdvdxdt) 



{edt + V ■ V^)(t)i = 7 ( e<5 



where 7(Z) = 27(Z) + Zi{Z), while Q,^ and Ql are defined in ([42]). 
Clearly, 7 has support in [0, 2], so that 



^l-l^^))(^^)^o(i 



On the other hand, the the fluctuation control (b) in Proposition 12.31 and 
the estimate (14. 6j) imply that 



7 I I : ) ) I \ ) = '->\^)L^{q^^{dx;L^{Mdv))) 
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Together with Lemma 14. 2^ these last two estimates lead to the following 
bound 

/ l-a/2\ 

{edt + V ■ = 0' ' 



\ / L^(Mdvdxdt) 

+ 0(l)L2(Li^^(dx;Li((l+kl)"'^/^^"«))) 



+ 



L]^^{dtdx;L'2({l+\v\)-f>Mdv)) 

Pick then a G (1, 2); the last estimate implies that ()4.7p holds for each 5 > 0, 
as announced. 



Stey 3. On the other hand, we already know from the fluctuation control 
(b) in Proposition 12.31 and (j4.5p that 

(4-8) 4)1 = 0(1)2,00 . 

Moreover 

(4.9) is locally uniformly integrable in the f-variable 

Indeed, for each ip G n L^(L^), one has 

2 



K 



< 1 1 { ^ ] \v\Mdxdv 



+ 



|v7|M(i^;(ix. 



The second term is 0(e"~^)||(/)||_Loo. Hence this term can be made smaller 
than any given 7] whenever e < eo(^)- Since e denotes an extracted subse- 
quence converging to 0, there remain only finitely many terms, say N = N{rj) 
that can also be made smaller that rj, this time by choosing ||^||l°°(li) 
smaller than c = c{N{rj),r]). As for the first term, it can be made less than 
T] whenever ||i;^||lc>o(2,i) < c'(?7), by Proposition 13. 2i Therefore 



K 



lipMdvdx 



< 2r] for each e and 5 > 



whenever ||</'||Lgo(Li) < min(c( A^(?7), r/), c'(r/)), which establishes ()4.9p . 

Applying Theorem lB.il (taken from [13j) in the Appendix below, we con- 
clude from (gS]), (gj]) and (gZl) that 
(4.10) 

for each 5 > 0, '/'f is locally uniformly integrable on R4. x x 



for the measure Mdvdxdt. 
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Step j. But, for each e, (5 G (0, 1), one has 

2 



Ve" + - 1 



1 — 7 eo 



2 



1 CI 



so that 



ln(5| / i,f{i\(Mdvdx)) 

by the fluctuation control (a) in Proposition 12. 31 This and (14.101) imply that 
(4.11) 

^/^^^+g;-i\' . „ .,^3..„3 



is also locally uniformly integrable on R+ X X 

for the measure Mdvdxdt. 

Because of the estimate (14. 5p in Step 1, we finally conclude that 

2 



(4.12) ( — — ] is locally uniformly integrable on R+ X R^ X R^ 



for the measure Mdvdxdt. 

Together with the control of large velocities in Proposition 13.21 the state- 
ment ()4.12p entails Proposition 14.11 □ 

Here is a first consequence of Proposition 14. H bearing on the relaxation 
to infinitesimal Maxwellians. 

Proposition 4.3. Under the assumptions of Theorem \2.4\ one has 



n-!— ^ >Q in L^^idtdx; L\{1 + \v\)PMdv) 



e e 
for each p < 2 as e ^ 0. 

Proof. By Proposition 14.11 the family 

2 



+ 

is uniformly integrable on [0, T] x K x R^ for the measure Mdvdxdt, for 
each T > and each compact K CH^. 
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On the other hand, ()3.ip and the the fluctuation control (b) in Proposition 
12.31 imply that 



fCT — 1 \ 
n^^— ^ j ^ in L\^JyMdvdxdl) 



and therefore in Mdvdxdi-measure locally on R_|_ x x R^. 
Therefore 

(1 + \v\f (^ ^"-^ - n ^"-^ ^ ^ in Ll^{dtdx-L\Mdv)) , 
which implies the convergence stated above. □ 



We conclude this section with the following variant of the classical veloc- 
ity averaging theorem [T2] , stated as Theorem IB. 21 in [13] . This result 
is needed in order to handle the nonlinear terms appearing in the hydrody- 
namic limit. 

Proposition 4.4. Under the assumptions of Theorem \2.4\ for each ^ G 
L'^{Mdv), each T > and each compact K C R^ 

f-T f 

\{ig,-f,){t,x + y) - {^g,j,){t,x)\'^dxdt^O 

JK 

as \y\ 0^, uniformly in e > 0. 
Proof. Observe that 

5.7. - 2^^^^^ = + - 2) ; 

since, up to extraction, 

+ 1)7e -2^0 a.e. and |(v^+ l)7e - 2| < 3 + \/2, 

it follows from Proposition 14.11 and Theorem lA.ll in the Appendix below, 
referred to as the Product Limit Theorem, that 

\/G~ — 1 

(4.13) ff,7, - 2-^— J > in Lf^^{dtdx; L^{Mdv)) 

as e ^ 0. 

This estimate, and step 1 in the proof of Proposition 14.11 (and especially 
the estimate (14. Sh there) shows that one can replace ge7e with +Gt- X. 
with Q > 1 in the equicontinuity statement of Proposition 14. 4[ 

Using (14. lip shows that, for each a £ (1, 2), the family 

is locally uniformly integrable on R+ x R^ x R^ 
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for the measure Mdvdxdt. In view of the estimate ()4.5p and Proposition 
13.21 we also control the contribution of large velocities in the above term, so 
that, for each T > and each compact K C R^, 

^ is uniformly integrable on [0, T] x K x 



for the measure Mdvdxdt. 

On the other hand, Lemma 14.21 shows that the family 



{edt + V ■ Vx) is bounded in L\^^{Mdvdxdt) . 

Applying then Theorem IB. 21 (taken from [13]) in the Appendix below 
shows that, for each T > and each compact K C R^, one has 



JK 



e / \ e 



dxdt 



as \y\ uniformly in e, which concludes the proof of Proposition 14.41 □ 



5. Vanishing of conservation defects 

Conservation defects appear in the renormalized form of the Boltzmann 
equation precisely because the natural symmetries of the collision integral 
are broken by the renormalization procedure. However, these conservation 
defects vanish in the hydrodynamic limit, as shown by the following 

Proposition 5.1. Under the same assumptions as in Theorem \2.4\ for each 
^ G spanjl, wi, t;2, ^3, |f P}, one has the following convergence for the con- 
servation defects J^eiO defined by : 

De(^) — > m L}^^{dtdx) as e ^ . 

Proof. For G span{ 1 , , f 2 , , bP}) the associated defect D£(^) is split as 
follows: 

(5.1) D,(e) = DKO+D?(0 

with 

and 

with the notation (j2T3|) and (pJ6]) . 
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That the term D^(^) vanishes for ^{v) = 0{\v\'^) as \v\ —>■ +00 is easily 
seen as follows: 



(5.2) 



r2 



< eO{K,)0{l) = 0(e|lne|), 



because of the entropy production estimate in Proposition 12.31 (d) and the 
choice oiK^ in ^M)- 



We further decompose D^(^) in the following manner: 

(5.3) D2(e) = D2i(0+Df(e) + Df(0 

with 



T V'-^e'-^el 



and, by symmetry in the v and v\ variables, 



Df (6 = -((e + 6)7e7ei7;7-i^ ^ : ^^^ TaoT 



The terms D^^(^) and D^^(^) are easily mastered by the following classical 
estimate on the tail of Gaussian distributions (see for instance [13j on p. 103 
for a proof). 

Lemma 5.2. Let GAr(z) he the centered, reduced Gaussian density in R'^, 
i.e. 



GNiz) 



(2^)^/2 



Then 



\z\PGN{z)dz ~ (2vr)-^/2|S^-i|i?H^-ie-^^ 



\z[^>R 



as R ^ +00. 
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Indeed, because of the upper bound on the colhsion cross-section in ()1.6|) . 
for each T > and each compact K C R^, 

l|DeH6llLi([0,T]x/f) 



< - 1 1 C 1 1 V 1 2 > if , Te vGe^^el 1 1 L 2 ( [0 , T] X ii: , L 2 ) 



el 



1/2 



<^||e^l|,p>^^(l + |z;|) 



/3||l/2 I 
L^(Mdv) I 



liai(i + i^il) 



/3||V2 



lLi{[0,T]xii',Li{Afidt)i)) 



f 2 



In the last right-hand side of the above chain of inequahties, one has obvi- 
ously 

From Young's inequality and the entropy bound (|2.19|) . we deduce that 



^,(1 + \v\f < (1 + Iwp) + 4 h{G, -l) + h 



1 + \v\ 



— O(l)Li([0,T]xi4',Li(Mdi;)) • 
Lemma 15.21 and the condition £^{v) = 0(|^;p) as \v\ — > +oo imply that 



v\r\\r\,.,.^. =0{K,^ e-2^^) = 0(ek/2|ine|V), 



on account of (12.430 . Thus 

(5.4) \\^fm\LHlo,T]>cK) = 0(6*^/'-'! Inel^) ^ 

for all .^(v) = 0(1^1^) as \v\ —>■ +oo as soon as k > 2. 

Next we handle D^^(^). Whenever ^ is a collision invariant (i.e. whenever 
^ belongs to the linear span of {1, wi, f2, fs, then C + ?i = + C'l; and 

using the {v,vi) — {v',v[) symmetry (j2.28p in the integral defining D^^(^) 
leads to 



Df(e) 



- y ^ + 6)7.7.i7e7.i — p 



Dfi(e)-Df^(e) 



232/ 



where 



Df Ho = h{ (c + ei)i|.|2+|,2|<^^7.7ei%%i 

and 



el) 



DP(e) = H + Cin\v\^+\v!\>Kje%lt%V 



/-/ iV ^'e^'el ~ V G^G, 



el 
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Then 



r2 



X ^||(^ + 6)l|^,|2+|^i|2</^,7e7€l%7£ll|L,-,,,,^,, 



= e-0(l).0(i^,)||7||ico 

so that 

(5.5) ||Df 1(011^1^ = O(eKe) ^ as e ^ 0. 

On the other hand, since G^ G [0, 2] whenever j{Ge) ^ 0, 

\\j^f\mL^^ < i6ii7iiiooi|i(^+ei)iH2+KP>KjiLi 



o 



+ bl^ + l^'lP)!! + 1^^ - ^^l|)'^l|,;|2 + |^i|2>i^JlLl(MMidt,dt,i) 
+ + \vi\'^)^~^'^^'^l\v\2+\vi\2>KML^MMidvdvi) 



0^]0 e 



0+6 



SO that 



(5.6) 



as e ^ , 



for k > 6, by a direct apphcation of Lemma 15.21 in x Rj^^ — i.e. with 
N = 6. 



Whereas the terms D^(^), D^^(^), and DP(^) are shown to vanish by 
means of only the entropy and entropy production bounds in Proposition 
2.31 (a)-(d) and Lemma 15.21 the term D^^(^) is much less elementary to 
handle. 

First, we split D^^(^) as 



Df (0 = -(c%ii-%i)VG:Gr, '^^'~'^^' 



+ -<|^C(7e7el(l - 7e) + 7e7el%(l " 7^l)) 

=Dfno+Dr(o. 



5 Y (-Te<-^el 
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For each T > and each compact K C R^, the first term satisfies 

v221. 
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<2Cb 



-{l-%i)VG7i{l + \vi\f'' 



\leyGl\\L° 



L'2{[0,T]xK;L^{R4idvi))) 



111^1(1 + l^'l)^/'llL^(Md.) 



el 



T2 



0(1) 



-(l-7.i)v^(l + l^i|)^/' 



L2([0,T]xA';L2(Mi(it)i))) 



provided that ^{v) = 0{\v\'^) for some m E N. 



Since supp(l - 7) C [|,+oo), then - JIJ2-1 



Gc ^ VV^'^ whenever 7^ 7^ 1, 



and one has 



1, 



Furthermore, as 

(5.7) |1 — 7e| < 1 + ||7'||l°° and 1 — 7^ — > a.e., 

the uniform integrabihty stated in Proposition 14.11 and the Product Limit 
Theorem (see Appendix A) imply that 



(5.8) 

Thus 
(5.9) 



^-11 



in L^([0,r] X K,L'^{M{l + \v\fdv)). 



Df (e)||Li([o,T]x/r)^Oase^O. 



Finahy, we consider the term D^^^(^): one has 

l|De^^(C)llLi([0,T]xif) 

2 / 

< - ( 11(1 - %)C\\l^[0,T]xK;LI) + 11(1 - 7ei)^lli2([0,T]xE-;L2) 



0(1) 
0(1) 



7.vaiiL,-^ 
1-7 



e 

1-7. 



(l + |z;|2 + |z;i|2) 



f 2 



LH{Q,T]xK;Ll) 



L'^{[0,T]xK;L'^{Mdv)) 

where the first equahty uses the {vvi) — {v'v'i) symmetry in ()2.28p . 
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Since supp(l — 7) C [|, +00), < ^ whenever % / 1, one has 



ii-7ep < V2 \i-%\vg:-i 



n- 



/GT-l 



By dSZI) and ([SSI) 

(5.10) E^<i±^and^ 



in Lf^^{dtdx,L^{Mdv)) 



since ^/Gl— 1 > \J?>l2 — 1 whenever 7^ 7^ 1, whereas by Proposition 12.31 (b) 
and Lemma |3. II 

^ = 0(l)ioo(j;^2(^9(JVMl)))) ) 



n- 



for ah g > +00. Then, 



|l-7. 



o(i; 



L}^Jdtdx,Li{Mdv)) 



for ah (7 < 2. In particular, for each r < +00, (^^{1 — + \v\Y) is 
uniformly bounded in Lf^^{dtdx, Lp'{Mdv)). By interpolation with ()5.10p we 
conclude that 

2 



(5.11) 



1-7. 



+ 



as e ^ 



L^{[0,T]xK:L'^{Mdv)) 

and consequently 

(5.12) Df ^(6 ^ in Ll^{dtdx) as e ^ . 

The convergences ([O]), dS^l), dSj]), (IET2D . (l53]) and ([52]) eventually 

□ 



imply Proposition 15. 1[ 

Remark. The same arguments leading to (j5.8p and to (jS.lip imply that, 
for each r G R, 

2 



(5.13) 



1-7. 



(1 + bl 



as e ^ . 



L2([0,r]xii';L2(Mdi))) 

6. Asymptotic behavior of the flux terms 
The purpose of the present section is to establish the following 
Proposition 6.1. Under the same assumptions as in Theorem\2.4[ one has 



F,(c)-2(c n 



/a - 1 



in Lj^^{dtdx) 
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as e — > 0, where C o-nd (" designate respectively either A and A or B and B 
defined by (KM) and (KM)- 

Proof. First, we decompose the flux term F^{() as follows: 



/a-i\ \ , 2/^ vg:-i 



= f1(c) + f2(c). 

We further split the term F^(0 as 



with 
(6.1) 



n 



+ n- 



^-1 



Fe''(C) = (C(i|.p</r.7.-i)(n 



Fi^(C)= C n 



/a - 1 



The term FP(C) is easily disposed of. Indeed, the definition (|2.32|) of 
the hydrodynamic projection 11 implies that ^n ^^^~^ -^ (1 + \v\Y is, for 
each p > 0, a (finite) linear combination of functions of v of order 0(|w|p+^) 
as \v\ +00, with coefficients that are quadratic in for ^ G 

{1, wi, ^2, ■^3, bp}- Together with Proposition l4.lt this implies that, for each 
r > and each compact K C R^, 

/ fo' _ 1 \ ^ 

(6.2) ( n-^— ^ j (1 + \v\Y is uniformly integrable on [0, T] x K x 

for the measure Mdvdxdt. On the other hand, 

l|i;|2<if,7e - 1 ^ and |l|^|2<j^^7e - 1| < 1 a.e. . 

Since C,{v) = 0(1^^) as \v\ +oo, this and the Product Limit Theorem 
imply that 

(6.3) Fl\C)^0 in LUdtdx). 
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The term F^"^((^) requires a slightly more involved treatment. We start 
with the following decomposition: for each T > and each compact K C R^, 



(6.4) 

l|Fe^(C)llLi([0,T]xA') < 



+ n-^— ^ — 



X 



e 

/a - 1 



n- 



e 



L2{[0,T]xii';L2{Mdi))) 



L2{[0,T]xK;L^{Mdv)) 



Since 7^ = 7(Ge) = whenever > 2, one has for each q < +00, 
2 

7£ 



In particular 



/a - 1 



L2([o,T]xii';L2(Mdt))) 



0(1) 



since C{v) = 0(|f as \v\ +00. This and (j6.2p imply that 



3.5) 



— ^ — + n-^^-^ — 



e e 
Using (16. 4p . ()6.5p and Proposition 14.31 shows that 
(6.6) Fl\C) ^ Ll,{dtdx) 

This and (j6.3p imply that 



Lf^Jdtdx;L2{Mdv)) 



0(1). 



(6.7) 



F.^(C)-(C n 



in Ljg^{dtdx) 



as e ^ 0. 



Next we handle the term F^((^). We first decompose it as follows: 

(6.8) +2(C^^^) + ^(C^^ 

= F2i(C)+Ff(C)+Ff(C). 
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Then, by ()2.20p and Lemma 15.2^ one has 
(6.9) 

21 2 
ll^e (C)IIl^{L2) < -||7l|L°°IICl|,;|2>A-J|L2(Mdi)) 



L^{L'^{Mdvdx)) 



On the other hand, for each T > and each compact K C R^, 
(6.10)^ 

\\^e {C}\\l^{[O.T]xK) 



< 2r^/^ 



7.-1 



L'^{[0,T]xK;L^{Mdv)) 



/a - 1 



Lf(L'^(Mdvdx)) 



as e ^ 



because of (|2.20j) and of (j5.13p . since C{v) = 0{\v\^) as |t;| +oo. 
Finahy, we transform F^^(C) as follows: 



Ff (C) = 2{C-^C 



2 C Q 



1 



Writing 



/a - 1 a/g: - 1 



+Q ( _ n^^^i^, + n- 



/a-i 



e e e 

and using the classical relation (see [2] for instance) 

Q((/>, (/>) = |/:((/>^) for each G Ker £ , 

we arrive at 
Q 



i/: I ( n- 



Thus 
(6.11) 



Ff (C) = ( C ( n 



+ 2(CQ 



/q - 1 

e 

x/a - 1 



|(CQ(\/G'„ 
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By continuity of Q (see ()2.35p ). 



< C'||C||L2{aAMv) 



e 

/a - 1 



Li([0,T]x_ft:) 



n- 



/a - 1 



— ^ — + n-^— ^ — 



L^{[0,T]xK;L2{{1+\v\)I^ Mdv)) 



L2([0,T]xi<';L2((l+|i)|)'3Mdv)) 

as e ^ 0, for each T > and each compact K C R^, because of (|6.2|) and 
Proposition 14.31 



Thus, by IKWh and (IHTTD 



(6.12) f2(c)-(c n 



/a - 1 



0. 



in Lj^^{dtdx) as e — > 0. 

The convergences (|6.7p and (|6.12p eventually imply Proposition 16. li □ 

7. Proof of Theorem 12.41 

Throughout this section U = U{x) designates an arbitrary compactly 
supported, C°^, divergence- free vector field on R^. Taking the inner product 
with U of both sides of (j2.48p gives 



(7 



dt J {vK.gele) -Udx- J F,{A) : V^Udx = j T>,{v) ■ Udx 

^ in Ll^idt) , 



by Proposition 15. 1[ Likewise, the energy equation ()2.49p and Proposition 
15.11 lead to 

5t(i(bP - 5)K,9e7e) + div, F,{B) = D,(i(|7;|2 - 5)) 

^ in Lj^^{dtdx) . 
By Proposition 16.11 one can decompose the fluxes as 



(7.2) 



(7.3) 
where 

(7.4) 



FM) = F™(A) + ^"(A) + o(l)ii^^(,i,.) 
F,{B) = F^o^^(i?) + F'^^iB) + o(l)z.i„^(,i,.) 



iiconv 



{A) = 2{A{n 



/a - 1 



FdifF(^) = -2(A^Q{y^„ 



/a - r ^ 
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while 

F™(5) = 2/5 (n- 

(7.5) \ V 

Fdiff(i?) = -2/s-2Q(\/G:,yG: 

Classical computations (that can be found for instance in [3]) using the fact 
that A is orthogonal in Lp'{Mdv) to Ker£ as well as to odd functions of v 
and functions of |f p show that 

pconv(^) = 2{A®A): (v ^'~^ 

In a similar way, B is orthogonal in L^{Mdv) to Ker C and to even functions 
of so that 

■pconv^ D\ _ 0/ D ^ D\ /„,v^^ ~ \ //I i„,i2 -na/G^— 1 



7.1. Convergence of the diffusion terms. The convergence of Fdiff(yl) 

and F^^ (i?) comes only from weak compactness results, and from the fol- 
lowing characterization of the weak limits. 



Proposition 7.1. Under the same assumptions as in Theorem \2.4\ one has, 
up to extraction of a subsequence 0, 



(7.6) 9e„^ 9, and ^ > q 

m w — L]^^{dtdx;L^{Mdv)) and in w — L?' {dtdxdij) respectively. 

Furthermore g € Lf^{L'^{dxMdv)) is an infinitesimal Maxwellian of the 
form 

(7.7) g{t,x,v) = u{t,x) ■ V + 6{t,x)^{\v\'^ — 5) , div^; n = , 
and q £ L'^{dtdxdfi) satisfies 

(7.8) / / qb{v - vi,uj)dojMidvi = ■ V^g = UA : V^n + B ■ VJ) . 



Proof. Proposition 12.31 (c) shows that 

{g^) is relatively compact \n w — L\^^{dtdx;L^{Mdv)) 
while (j2.2ip implies that 



is relatively compact in tf; — L^^dtdxd^) 



c2 



Pick then any sequence e„ ^ such that 



gen ^ 9, and 2 ' ^ 
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m W — L}^^{dtdx;L\Mdv)) and in w — L'^{dtdxdn) respectively. 

Step 1 : from (|2.22|) we deduce that 

^ginw- Ll^{dt,L\dxMdv)) . 

In particular, by Proposition 14.3^ g is an infinitesimal Maxwellian, i.e. of 
the form 

g{t, X, v) = p{t, x) + u{t, x) ■ V + 6{t, x)^{\v\'^ — 3) . 
Taking limits in the local conservation of mass leads then to 

divx{vg) = 0, 

or in other words 

divx u = 

which is the incompressibility constraint. 

Multiplying the approximate momentum equation (j2.48p by e 

edt{vK,gele) + ediv^F,(^) + \V x{\\v\\j,-t,) = eTi,{v) , 

using Propositions 15. l l and 16. Il to control 'D^{v) and the remainder term in 

and estimating F^°^^{A) and F^^^(74) by the entropy and entropy produc- 
tion bounds (l2:2nD - (l2:nD 

(^(n^^^^) ) = 0(i)mLr(4), 

we also obtain 

yx{\v\^g) = o 

or equivalently, since {\v\^g) = 3{p + 6) £ L°°(R+; L2(r3)), 

P + e = o, 

which is the Boussinesq relation. One therefore has (j7.7p . 
Step 2 : start from (14. ip in the proof of Lemma 14.21 : 

Recall that 

(7.9) Q\^^inL^[Mdvdxdt) . 
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Next observe that 

-,2 \fG't 
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—— \J G' G'^ — ^ GfGfX 
G^i— — b{v - vi,uj)dLoMidvi 



V^^^TG. 

Proposition 14.11 implies that 

^/G', 1 in Lf^^{dtdx; ^^((l + \v\fMdv)) as e ^ ; 
this and the second Hmit in (j7.6p imply that 



1 5 b{v — vi,iJ)dujMidvi 



qb{v — vi,uj)dujMidvi 



in w — Llg^{dtdx; L}{Mdv)) as n ^ +oq. Since on the other hand 



, . 1 a.e. as e — > with < , 

+ G, - Ve" + G, 

we conclude from the Product Limit Theorem that 

-i2 



(7.10) 



Qi 



qb{v — vi,uj)dujMidv 



in w — Ll^^{dtdx; L^{Mdv)) as n — > +oo. 
By (gSD, (imH and TO 



in w - L}^^{dtdx;L^ {{I + \v\'^)Mdv)) whenever a e]l, 2 [. Using ([Tj]), (1710]) 
and the convergence above, and passing to the limit in (j4.ip as ^ leads 
to 

b{v — vi,u!)dLoMidvi = • Vxd , 



which is precisely the first equality in (|7.8p . Finally, replacing g by its 
expression ()7.7p in the formula above leads to the second equality in (|7.8p . 

□ 

Since A and B G L'^{aMdv), the second limit in (|7.6p and identity (|7.8p 
show that 



Fi}^iA) = -2(A^QUG,,VG 



1 



(7.11) 



-{Ai^ A) : VxU = -iy{VxU + {VxU 



Fdiff(i3) = -2(i?^Q(Va,VG, 



1 



-{B(g)B) : = -kVx 



in w — LF'{dtdx) as e — > 0, because of the divergence- free condition in (|7.7p . 
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7.2. Convergence of the convection terms. The goal of this section is 
to establish that 



u (X" n : VxU dx 



and div^F™(S) ^ |div^.(ue) 

in the sense of distributions on R!^ and on R;^ x R^ respectively. 

First, we replace Y^'^''^^ [A) and F^'-'^^^(S) by asymptotically equivalent 
expressions. 

Indeed, because of (j4.13p 

{vQele) - 2(v ^'~^ '^ ^ in Ll^{dtdx) 

and 

- 1)5,7,) - 2{{\\v? - 1)^^^^^ - in Ll,{dtdx) . 

On the other hand, (7,7, is bounded in Lf[L'^{Mdvdx)) while f 1|^|2>^^ 
and — l)l|^|2>j^^ ^ in L'^{Mdv); therefore 

in Lfg^{dtdx). Therefore 
(7.12) 

= {vK,9ele)®^ - \\{vK,gele)\^I + o(l)il^^(did^) , 

while 

F™(i?) = {B(E)B)- {vK,ga,){{\\v\^ - l)K,9ele) + 0(1) rl ^atd.) 

(7.13) 

Furthermore, since g^^ g weakly in Lj^^{dtdx; L^{{1 + \v\'^)Mdv)) while 

VKtJe V and (^1^1^ — l)_R'e7e (ibP ~ 1) ^-^-j and 
\vkM + - 1)kM < Cil + \v\^) 

one has by the Product Limit Theorem 

^ ■ ^ ((iH^-l),,,„7e„5eJ-((ibP-l)5)=^ 

in It; — L\^^{dtdx). In fact, these limits also hold in ti; — Lf^^{dtdx) since the 
family (7, 7, is bounded in L^{L'^{Mdvdx)). 

Taking limits in ()7.12p and (j7.13p . which are quadratic functions of the 
moments, requires then to establish some strong compactness on [{CK^Qele))- 
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7.2.1. Strong compactness in the x-variable. Applying Proposition 14.41 with 
^ = v and = P — 5) shows that, for each T > and each compact 

K c R^ 

r [ - 5)g,„j,Jit,x + y) - {^i\v\^ - 5)g,„^,Jit,x)\^dxdt 

Jo J K 

+ [ [ \{vg,r.l^J{t,x + y) - {vge„j,J{t,x)\'^dxdt ^ 
Jo Jk 

as \y\ — > uniformly in n. An easy consequence of the above convergence 

properties is that 

(7.15) 

/ / \{\{\v? - 5)K.^ger.7e„)it,x + y) - {^i\v\^ - 5)K.^g,„j,Jit,x)\^dxdt 
Jo Jk 

+ / \{vK,„ge^leJ{t,X + y)- {vK,„ge^leJ{t,x)\'^dxdt^O 

Jo Jk 
as |y| — > uniformly in n. 

In order to study the convergence of F^{A), we need some similar state- 
ments for the solenoidal and gradient parts of {vKe„9enlen) , since the first 
one is expected to converge strongly in Lf^J^dtdx). 

The difficulty comes then from the fact that the Leray projection is a 
non local pseudodifferential operator, in particular it is not continuous on 
LlMx). 

Introducing some convenient truncation x in x and using the properties 
of the commutator [x, P], one can nevertheless prove the following equicon- 
tinuity statement (see Lemma IC.ip : for each compact K C R'^ and each 
r > 0, one has 

(7.16) / / |P(t;i^,„(7.„7.J(t,a; + y)-P(fif,„5.„7eJ(t,^)|'dxdt^0 
Jo Jk 

as \y\ — > 0, uniformly in n. 

7.2.2. Strong compactness in the t-variable. As we shall see below , the tem- 
perature fluctuation (|(|vp — 5)_ft'^5e„7e„) and the solenoidal part P^VK^de'^e) 
of {vK^Qtlt) are strongly compact in the t-variable. However the orthogonal 
complement of P{vKe9tlt) — which is a gradient field — is not in general 
because of high frequency oscillations in t. 

Proposition 7.2. Under the assumptions of Theorem \2.4\ one has 

P{vK,„ge„7e„) {vg) = U 

(i(|r;|2-5)K.„<7.„7.J-(^(H'-5)<7) = i^ 
in C(R+; w — L^) and in Lf^^{dtdx) as n ^ +oo. 
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Proof. The conservation law ()7.ip implies that 

(7.17) dt [ {vK,^g,„j,„) ■ Udx = 0(1) in LUdt) 

for each compactly supported, solenoidal vector field U G C°°(R^), since we 
know from Proposition 16.11 together with the bounds ()2.2ip and ()2.20p that 
Fe„(^) is bounded in L\^^{dtdx). 

In the same way, the conservation law (17. 2p implies that 



(7.18) dt{\{\v\'' - 5)k,„5.„7.J = 0(1) in Ll^{dt- W-l^\n')) . 

Also, we recall that g^'^^ is bounded in B{^^]L?'{Mdvdx)) — where 
B{X,Y) denotes the class of bounded maps from X to y — because of 
the entropy bound ()2.20p . Indeed, since 7e = whenever > 2, one has 



(7.19) b.„7ej < la<2^^ < (1 + V2) '^ 



In particular, one has 
(7.20) 



{vK.,.9e„7eJ = 0(1) mB{R+;Li,, 



iUM' - 5)i^,„5e„7.J = 0(1) in B{R+;L 



X) 1 



Since the class of C°°, compactly supported solenoidal vector fields is 
dense in that of all 1? solenoidal vector fields (see Appendix A of [19j). 
(I7:20]l and (ITTfl) imply that 

(7.21) P{'^Ke^genlen) is relatively compact in C(R+;it; — -L^(R^)) , 

by a variant of Ascoli's theorem that can be found in Appendix C of [19]. 
The same argument shows that 

(7.22) {\{\v\^ -^)Ke^ge„ler,) is also relatively compact in C(R+;t(; - lI) . 

As for the Lf^J^dtdx) compactness, notice (|7.2ip - (|7.22p imply that 
P{'VK,„ge„le„) -*^X5 is relatively compact in Lf^^{dtdx) 
(i(|up - 5)ii-^^fi(,„7,„) is relatively compact in Lf^^{dtdx) 

where X5 designates any mollifying sequence and * is the convolution in the 
j;-variable only. Hence 

P{vK,„ge„le„) ■ P{vK,^ge„leJ -k XS ^ Pu ■ Pu -k XS 

{\{\v\^ - 5)i^.„5e„7.J(^(bl' - 5)k.„5.„7.J - (1^) [¥^X5) 
m. w — L\^^{dtdx) as n ^ oo. By ()7.15p - ()7.16p . 

P{vK,„genleJ *X5 ^ P{vK,^ge„leJ 

(7.23) 
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in Lf^^{dtdx) uniformly in n as (5 ^ 0. With this, we conclude that 

\{\{\v? - 5)/r.„5.„7.JP - (^^) in «^ - Ll,{dtdx) 
which implies the expected strong compactness in Lf^J^dtdx). □ 

7.2.3. Passing to the limit in the convection terms. As explained above, 
P{i'K^„9e„j€n) is strongly relatively compact in Lf^^{dtdx); however, the 
term {vK^^Qenltn) itself may not be strongly relatively compact in Lf^^{dtdx) 
— at least in general. For that reason, on account of (j7.12p . it is not clear 
that 

pConv(^)^^^^_ i|^|2^_ 

Likewise P — 5)Ke„ge„'ye„) is strongly relatively compact in Lf^^{dtdx), 
and, on account of (I7.13p . it is not clear that 

pconv(^) ^ 5^0 

as one would expect. 

What we shall prove in this section is 

Proposition 7.3. Under the assumptions of Theorem \2.4\ one has 

[ V^U : Ff'^™(^)dx ^ / V^U ■.u(S)udx 

in the sense of distributions on for each solenoidal vector field U £ 
C^°^(R3;R3), and 

div,.Ff«™(S)^|div,(ne) 
in the sense of distributions on R^ x R'^. 

The proof of this result relies on a compensated compactness argument 
due to P.-L. Lions and N. Masmoudi [21] and recalled in Appendix A (The- 
orem A. 2), and on the following observation: 

Lemma 7.4. Let 6 > 0, and ^ G C^(R'^) be a bump function such that 
supp^ C 5(0, 1) , ^>0, and j^dx = l; 

let ^s{x) = S^^xi^/^) '^''T'd ^5 = ^5*^5'*^Cs- Denote by Q = I — P the orthog- 
onal projection on gradient fields in L^(R'^;R^). Under the assumptions of 
Theorem \2.4l ''^^ has 

edtQ{\s * {vK,„ge„-feJ) + V^Xs ★ {^\v\K^^ge„%J ^ 

edtXs -k {l\v\K,^ge„'ye„) + | div^ Q{Xs -k {vK,„ge„7e„)) 

in L]^^{dt;Hl^{H?)) for each s > 0. 
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Proof. The second convergence statement above is obvious: indeed, consid- 
ering the truncated, renormahzed energy equation (I2.45P with ^ = 
and applying the mohifier leads to 

edtXs -k {l\v\K^^^ge„%„) + f div^. ^(A^ -k (uK,„ffe„7.„)) 
= -fediv^ Xs*F,{B) + leXs * D,(|t;|2) . 

It follows from Proposition I6.H the entropy bound (|2.20p and the entropy 
production estimate ()2.2ip that Ff:(B) is bounded in Lj^^(dtdx); this and 
Proposition 15.11 eventuallv entail that the right-hand side of the above equal- 
ity vanishes in Ll^^{dt; Hf^^(R^)) . 

The first convergence statement above is much trickier. Start from the 
analogous truncated, renormahzed momentum equation ()2.45p with ^ = v: 

(7.24) edt{vK^^ge^7eJ+Vxi{\v\K,^9e„leJ = "6 div, F, (A) + eD, (t;) 

Applying Q to both sides of the equality above is not obvious, because we 
only know that the right hand side vanishes in L\^^{dt] Wi^^'^^R?)), while Q 
is known to be continuous on global Sobolev spaces only. 

However, Q = V^A"^ div^; is a singular integral operator whose integral 
kernel decays at infinity. More precisely, we shall use Lemma IC.21 together 
with the following estimates on the right hand side of (j7.24p : 

Lemma 7.5. One has 

eFe(A) and eT)^{v) ^ in L]^^{dtdx) 
as e ^ 0. Furthermore, 

eF,{A)=0{l)L^(^Ll). 
eD,{v) = 0{e^Ky')Li^ + 0(V^)L?(Li) + ^(1)^?,. • 

Note that the convergence statement in Lemma 17.51 is a simple conse- 
quence of Propositions 15.11 and 16.11 (already used in the derivation of the 
Boussinesq relation in paragraph I7.ip . 

Then let us postpone the proof of the global estimates, which is based on 
the entropy and entropy production bounds (j2.20p - (j2.2ip . and conclude the 
proof of Lemma 17.41 

Define Q = (,5* S,5- One has then 

edtQiCs* {vK.gele)) +V^C5* {^\v\lje%) = -QiCs * {^F .{A) * ^ ^s)) 

+Q{C5*{eF>,{v))). 

For each 5 > fixed, 

Qi^s * (.eF,{A) * V^s)) ^ in L}^^{dtdx) as e ^ 

by the first convergence result in Lemma 17.51 and Lemma |C.2[ 
Next decompose 

eD,(7;) = DO(t-)+D;(z;) 
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with 

B'Av) = 0(1)l?,, and D;(^;) = Oie^K}/^)^!^ + . 
Thus, one has 

Cs^-D'M^O in Ll,{dt;Ll) 

so that 

Q{Cs*-D[{v))^0 in Ll,{dt;Ll) 

as e ^ 0, by the L^-continuity of pseudo-differential operators of order 0. 
Finahy, since D^(f) — > in Lj^^{dtdx) and is bounded in L^^, it follows 
from Lemma IC.21 that 

QiCs * DO(w)) ^ in Ll^idtdx) . 

Eventually, we have proved that 

edtQiCs * {vK.Qele)) + * {\\v\'K,9ele) ^ 

in L]^^{dtdx) as e — > 0. Therefore, denoting A5 = ^5 * ★ ^s, one has 

edtQ{\s * {vK.gele)) + V^A^ -k {^\v\Kje7e) ^ 

in Lj^^idt; H[^^(R^)) for each s > as e ^ 0. □ 
Let us now turn to the 



Proof of LemmaYT^ First, 51,7^ = 0(1) in Lf{L'^{dxMdv)), while A E 
L'^{Mdv): hence 

Next decompose eY)^{v) as 

eD,(t;) = ri + r2 + r3 

where 



T2 = 2(( vkA, 



n = 2^VKneVG'e [VgTi " i) ^ (^^^ - v^g^oT) 



Since ^ — \/G^G^j is bounded in (see (I2.2ip ) , one has 

Ti = 0(e2i^V2)^^,^^ . 

Likewise, 7e\/G7 = 0(1) in Lf^^ ^ and u € L'^{dfi), so that 

r2 = o(ik. . 
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The same argument is used for T3, except that one has to control the 
terms v [y/Gei — l) instead of v in L^^. By Young's inequality, 



{l + \vi\ 



G.i-l <{l + \vi\)\G 



el 



< - (/i(Gei - 1) + r (e(l + |^;i|))) 

< -hieg,i)+eh*{l + \vi\) 
e 



The 3rd inequality above comes from the superquadratic nature of h* . In- 
deed 



£ — / n 



n>2 



SO that 

h* (Ap) < A^/i* (p) , for each p > and A G [0, 1] . 
With the upper bound on J h{v — vi,uj)duj, this shows that 



l^sl < \\VK, 



\L^{{l+\v\)0Mdv) 



lleGeWlSf \\V Gel — 1| 



L2{{l+\vi\)I^Midvi) 



~^ G'eG\\ — \l G^Gel^ 

Combining the previous results leads to the expected estimate for T)^{v). □ 



At this point, we conclude this section with the 



Proof of Proposition 1.3. First, we apply the compensated compactness ar- 
gument for the acoustic system in |21| — see also Theorem lA.2l — to conclude 
from the statement in Lemma 17.41 that 

V,,C/ : Q{\, * {vK,^g,^leJ)''^dx ^ 

div^(A5 ★ {l\v\K,^ge„1e„)Q{Xs * {vK,^ge„7e„))) ^ 

in the sense of distributions on and x R^ respectively, for each 
divergence- free vector field U G C^(R^;R'^). 

On the other hand, the compactness property in the x-variable stated in 
Proposition 14.41 and (j7.23p implies that 

Qih * {vK,„ge„leJ - Q{{vK,„ge„leJ ^ 
* {l\v\K^^gen7eJ - {l\v\K^^gen7eJ ^ 
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in Lf^^{dtdx) as 5 ^ 0, uniformly in n. Therefore, one has 



(7.25) 







in the sense of distributions on H*^ and H*^ x R'^ respectively, for each 
divergence- free vector field U £ C^(R^;R'^). 

Also, we recall from Proposition 17.21 and ()7.14p that 

P{vK,^9e,,leJ) u strougly in Lf^^{dtdx) , 
Q{vK,„ge„l,J weakly in Lj^^{dtdx) . 

Therefore, for each compactly supported, C°° solenoidal vector field U , one 
has 



R3 



+ 



R3 



V.U:{P{vK,^g,^le„)r^dx 
V.,U:{Q{vK,Mejr'dx 



+ 



R3 



V^C/ : {P{vK,„ge,,7eJ <^ Q{vK,„ge„len))dx 



+ / V^C/ : iQ{vK,„ge„7eJ <^ P{vK,Qe^le^))dx 



R3 



/ 

JR3 



VxU : u (8> udx 



in the sense of distributions on R!j.. Together with (j7.12p . this implies the 
first limit in Proposition 17. 3[ 

On the other hand. Proposition 17.21 and (|7.14p imply that 

- l)/f.„5'e„7€„> G strongly in Lf^^{dtdx) , 
{\v\K,ge„7e„) weakly in Lf^^{dtdx) . 



Hence 



divx{{{h\v\ 



,9e„7e„){vK,„ge„7e 



= div^(((ibl^ - ^)K,^genleJP{vK,„ge„le„ 

+ ^ div^. {{\v\K^^ge„%JQ{vK,„ ge„ 7e„ 

+ TE diVx ( ( b I ge„ 7e„ ) P{VK,„ ge„ 7e„ 
+ div^(((i|v|2 - '^)K,„gtn7eJQ{vK,„ge„%„ 

— > diVxiuO 

in the sense of distributions on R^ x R^. With (j7.13p . this entails the second 
statement in Proposition 17.31 □ 
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7.3. End of the proof of Theorem 12. 4L At this point we return to the 
renormalized, truncated momentum and energy conservations in the form 
(EH) and ra. 



Asymptotic conservation of momentum : by using the convergence prop- 
erties in (jT.lip and Proposition 17.31 with the decomposition (j7.3p . one sees 
that, for each C°°, compactly supported, solenoidal vector field U, 

/ V^U :F^^{A)dx ^ / V : u®udx-u \ V : {V ^u + {V xuf)dx 
Jr3 Jr3 JR3 

in the sense of distributions on R^, while 

div^F.jS) ^ div^(ne) - k/\J 

in the sense of distributions in x R'^. Furthermore, since div^; u = 0, one 
has 



/ V^U : {V^ufdx = I V^.(div^. U) ■ udx = 
Jr3 Jr3 

ich solenoidal test vector field U, so that 

/ VxU : F^^{A)dx / V^U : u (g) udx - u I V^U : V^udx 

JR3 Jr3 Jr3 



in the sense of distributions on Yi*_^. 

On the other hand, by Proposition 17.21 



U ■ {vK,„ge„leJdx ^ U ■ udx 

R3 JrS 

uniformly on [0, T] for each T > 0. In particular, for t = 0, one has 

[/ • .dx = lim / U-P\- I vFTdv\dx= I U-u'^'dx. 



R3 '* ° ^^Ojj^a J-R? ' J Jr3 

Therefore, u satisfies 

dt i U ■ udx - / VxU : u ® udx + v j V^U : V^udx = , t> 0, 
Jr3 Jr3 Jr3 

Asymptotic conservation of energy : likewise, 
in C(R+;t(; — L^). In particular, for t = 0, one has 



= w-lhn - [ {l\v\^-l)Frdv 



Therefore, 9 satisfies 

dtO + diY^iuO) - kA^O = , xeK^ , t>0, 



\t=o 
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Notice that one has also 

1 r 

vFe^dv = {vgej u 
■1 - - M)dv = {{\\v\'' - l)ff,J ^ 9 

weakly in Lj^^{dtdx), because of (|7.6p and (|7.7p . 

Asymptotic energy inequality : by Proposition 17.11 and (j2.22p . one has 



1 



^ginw- Ll^{dt,L\dxMdv)) 



and 

1 



2 - \/Ge„iaj - g in iz; - L^dtdxdfi) . 

Then, by (l2:20]l and ([MI]), 



// 



2 



Mg\t,x,v)dxdv < ]\m A I I M ( ^^^^ — 3^ ) {t,x,v)dvdx 



< hm ^i/(F,jM)(t), 

and 

2 

^ ri r r I \ I ^ Gi — \/ Gf„ i Gf„ \ 



q dsdxdfj, < lim / / / 2 I dsdxdfj, 

n^oo J J J \ 

< lim ^ / / £{FeJdsdx 
Explicit computations based on the limiting forms (j7.7p and (j7.8p of 
g and jj qb{v — vi,uj)dujMidvi 
and using the symmetries of q under the d/i-symmetries imply that 

//M/(,.,.M..„ = /,|„|^(t,.) + ||.|^(,.)M., 

while 

(see Lemma 4.7 in [3j for a detailed proof of these statements.) 
Taking limits in the scaled entropy inequality 

\H{F,\M){t) + \ f [ £{F,){s,x)dxds < \H{Fr\M) 
e ^ Jo J £ 
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entails the expected energy inequality 

/ {hu{t,x)\'^ + l\e{t,x)f)dx + [ I {v\V^u\'^ + ^K\VJ\^)dxds 

J-R? Jo Jr3 

<ljmirH(FnM) 
With this last observation, the proof of Theorem 12.41 is complete. 



Appendix A. Some results about the limits of products 

For the sake of completeness, we recall here without proof some classical 
results used in the present paper to pass to the limit in nonlinear terms. 

The first one is due to DiPerna and Lions [8], and is referred to as the 
Product Limit Theorem in [3j : 

Theorem A.l. Let ^ he a finite, positive Borel measure on a Borel subset X 
ofR.^ . Consider two sequences of real-valued measurable functions defined 
on X denoted ipn and ipn- 

Assume that {ipn) is bounded in L°^{dfi) and such that ipn ^ o-e. on 
X while ipn ^ in w — L^{dfi). Then the product 

fn'4'n ~^ fi^ in L^{diJ,) weak. 



The second one is due to Lions and Masmoudi |21j . and can be viewed as 
a compensated compactness result. It states that (fast oscillating) acoustic 
waves do not contribute to the macroscopic dynamics in the incompressible 
limit : 

Theorem A. 2. Let 0. Consider two families (ipe) and (V^V^e) bounded 
in L'^^{dt, Lf^^{dx)) , such that 

1 1 

dtife + -^xlpe = -Fe, 

e e 
c2 1 

€ € 

for some F^,G^ converging to in Lj^^{dt, Lf^^{dx)). 
Then the quadratic quantities 

PVx ■ ((V,Ve)®') and V, • (p.Vx^Pe) 

converge to in the sense of distributions on H*^ x R"^. 
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Appendix B. Some regularity results for the free transport 

operator 

The main new idea in our previous work on the Navier-Stokes hmit of 
the Boltzmann equation [iSj was to improve integrabihty and regularity 
with respect to the x variables using the integrabihty with respect to the v 
variables. 

This property is obtained by combining the velocity averaging lemma 
[111 I12j with dispersive properties of the free transport operator fB]. 

We state here two results of this kind used in the present paper, whose 
proof can be found in [H] or [13] . 

The first such result, based on the dispersive properties of free transport, 
explains how the streaming operator transfers uniform integrabihty from the 
V variables to the x variables. 

Theorem B.l. Consider a bounded family [ipf,) of L'^^{dt,L\^^{dxdv)) such 
that {edt + v • VxOV'e is bounded in Lj^^{dtdxdv) . Assume that [ipf,] is locally 
uniformly integrable in the v-variable — see Proposition \3. ^1 for a definition 
of this notion. Then {il^^) is locally uniformly integrable (in all variables t, 
X and V ). 

The second one, which is based on the classical velocity averaging theorem 
in [TTl [12], explains how this additional integrabihty is used to prove a 
averaging lemma. 

Theorem B.2. Consider a bounded family [ip^) of Lf^^{dtdx, L'^{Mdv)) 
such that {edt + v • Vx)^e is weakly relatively compact in Lj^^(dtdxMdv). 
Assume that (Iv'eP) is locally uniformly integrable with respect to the mea- 
sure dtdxMdv. 

Then, for each function ^ = ^{v) in L'^{Mdv), each T > and each 
compact K C R^, 



^■5 \y\ —^0 uniformly in e. 



Appendix C. Some regularity results for the Leray projection 

One annoying difficulty in handling incompressible or weakly compressible 
models is the nonlocal nature of the Leray projection P — defined on the 
space L^(R'^; R'^) of square integrable vector fields, on the closed subspace of 
divergence- free vector fields. By definition, P is continuous on L^(R'^;R^), 
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as well as on H^(Il^; R^) — since P is a classical pseudo-differential operator 
of order 0. However, P is not continuous on local spaces of the type L^^^{dx). 
Here is how we make up for this lack of continuity. 

A first observation leads to a local L^-equicontinuity statement provided 
that some global bound is known to hold. 

Lemma C.l. Consider a sequence of vector fields {tpn) uniformly bounded in 
Lf{Lp'{dx)). Assume that, for each T > and and each compact K C R^, 

/ / \il)n{t,x + y)-il)n{t,x)\^dxdt^Q as\y\^Q 
Jo Jk 

uniformly in n. 

Then, for each T > and and each compact X C R^, 

/ / \PMt,x + y)-PMt,x)\'^dxdt^O as\y\^0 
Jo Jk 

uniformly in n. 

Proof. For each 6 G (0, 1) and i? > 0, let x = x(x) be a C°° bump function 
satisfying 

x{x) = 1 for \x\ < R, x{x) = for |x| > i? + (5 , 

0<X<1, |x'|<2/<5. 
Obviously, for \y\ < 1, one has 

\x{x + y)^pn{t, x + y) - x{x)ipn{t, x)\^dxdt 



JR3 

T 



<2 [ [ x{x + yf\M-t,x + y)-Mt,x)fdxdt 

Jo iR3 

+2 [ [ \x{x + y)-xix)\''\Mt,x)\''dxdt 
Jo J-R^ 

<2 I [ \i;nit,x + y)-Mt,x)\''dxdt2(^] \y\''T\\i,n\\Lr^Ll) 

Jo J\x\<R+2 \0j ' ^ 



l\x\ 

SO that 

|2 



\x{x + y)V'n(i, x + y) - x(a;)V'n(i, x)\ dxdt 

Jn^ 

as |y| ^0 uniformly in n, since ipn is bounded in Lf°{L'^{Mdvdx)). 
Consider next the decomposition 

xP = Px + [x,^] 

where x denotes the pointwise multiplication by the function X) which is 
another pseudo-differential operator of order on R^. In particular, [x^P] 
is a classical pseudo-differential operator of order — 1 on R^. 
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With this decomposition, we consider the expression 

r-T f 

Xix + y)Pipn{t,x + y) - x{x)Pipn{t,x)\'^dxdt 

J\x\<R 

<2 [ [ \P{x^n){t,x + y)-P{xiJn){t,x)\^dxdt 

Jo J\x\<R 

+2 I [ \[x,P]Mt,x + y)-[x,P]Mt,x)\^dxdt 

Jo J\x\<R 

Because P is an L-^(dx)-orthogonal projection, the first integral on the right- 
hand side of the inequahty above satisfies 

T r 

/ \P{x'4'n){t,X + y) - P{x^n){t,x)\^dxdt 

J\x\<R 

< f [ \P{x^n){t,X + y) - P{xiJn){t,x)\''dxdt 

Jo Jr^ 

T 

2, 



< / \x{x + y)'4^n{t,x + y)-x{x)'(pn{t,x)\ dxdt^O 
Jo Jr3 

as \y\ 0, uniformly in n. On the other hand, because [x, -P] is a classical 
pseudo-differential operator of order —1 on (see [25], §7.16, on p. 268): 
therefore [Xt P] maps L^(R^) continuously into H^(R.'^). This implies in 
particular that [x,P]ilJn is bounded in L°°(Il-^-■, H^CR?)) so that, for each 
R>0, 

T f 

/ \[x,P]^IJn{t,x + y) - [x,P]Mt,x)\^dxdt ^ 

J\x\<R 

as |y| — > 0, uniformly in n. Hence 
rT r 

\x{x +y)PTpn{t,x +y) - x{x)P^n{t,x)\'^dxdt 

J\x\<R 

as \y\ 0, uniformly in n. Assuming that R > 2, that the parameter 6 in 
the definition of x satisfies 6 E (0, 1) and that \y\ < 1, we conclude that 

rT r- 

/ \Pi)n(t,X + y)-Pi)n{t,x)\^dxdt^{) 

>0 J\x\<R-2 

as \y\ — > 0, uniformly in n, for each R > sufficiently large. □ 



A second observation provides continuity of P in Lj^^ under some appro- 
priate global bounds. 

Lemma C.2. Let tp^: = tp^{t,x) £ R'^ be a family of vector fields such that 
V'e — > m L\^^{dtdx) and t/j^ = 0(1) in Lj^^{dt; L'^) . Let be a mollifying 
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sequence on R'^ defined by ^s{x) = 6^'''^{x/5) where ^ G C^(R'^) is a hump 
function such that 

supp«CB(0,l). «>0,«„./cdx = l. 

Then, for each 5 > 0, 

Q{^5 ★ Ve) ^ in Ljg^(dtdx) as e — > . 
Proo/. Let x G ^^(R^). Then 

nxix)\Qi^S*A)it,x)\dxdt = [ [ x{x)^{t,x)-Q{$,s-kil)e){t,x)dxdt 
where is any measurable unit vector field such that 

^{t.x) = ^^^^!l^(t,x) whenever Qi^s * ^e)it,x) 7^ 0. 
Hence ^ 

Jo Jn? 

= - / divj; (x^) diy^a {Cs * A) {t, x)dxdt 

Jo 7r3 

Let Q{x) = -^j^^s be the convolution kernel corresponding to — Va;A~^; for 
R>0, denote Gr{x) = ^(x)l|j.|<ij and g^{x) = ^(x)l|j.|>ij. Thus 

/ / x{x)\Qi^S*i'e)it,x)\dxdt 
Jo 7r3 

= f f gR*{x^){v^s)*Mt,x)dxdt 

Jo Jn^ 

+ I I g^*{x^){V^s)*Mt,x)dxdt 
Jo Jr3 

We have used here the following symplifying notation: if a and b are two 
vector fields on R^, we denote 



a-kb{x)=l a{x - y) ■ b{y)dy 
Jti» 



where • is the canonical inner product on R^. 

Observe that = 0{1/Vb,) in L^, while x^ G Lf{Ll.) (since \n\ = 1 
and supp(x) is compact). Hence 

g« * ixn) = 0{l/^) in Ll^{dt- Ll) 

and (V^^) * V'e = 0(1) in L]^^{dt; L^) for each (5 > since ■^g = 0(1) in 
Lj^^{dt; L"^) . Hence the second integral is 0(l/\/P) for each 6 > 0. 

Next Gr = 0{R) in and thus Qr (x^^) = 0{R) in since \n\ = 1; 
moreover, 

supp^ {Gr * (x^^)) C supp(x) + S(0, i?) 
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is bounded for each R > 0. On the other hand V^^ :A- ^ in Lj^^{dtdx), 
so that the first integral vanishes as e — >^ for each S > and each R > 0. 
Passing to the hmsup as e — 0"^, then letting i? — 0"^ leads to the announced 
result. □ 
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